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In the preceding paper, a rigorous three-dimensional relativistic equation for two-gluon bound 
states was derived from the QCD with massive gluons and represented in the angular momentum 

■ representation. In order to apply this equation to calculate the glueball spectrum, in this paper, 
' the equation is recast in an equivalent three-dimensional relativistic equation satisfied by the two- 
gluon positive energy state amplitude. The interaction Hamiltonian in the equation is exactly 
derived and expressed as a perturbative series. The first term in the series describes the one-gluon 

^ I exchange interaction which includes fully the retardation effect in it. This term plus the linear 

^ ■ confining potential are chosen to be the interaction Hamiltonian and employed in the practical 

calculation. With the integrals containing three and four spherical Bessel functions in the QCD 
CO , vertices being analytically calculated, the interaction Hamiltonian is given an explicit expression in 

the angular momentum representation. Numerically solving the relativistic equation with taking 
the contributions arising from the retardation effect and the longitudinal mode of gluon fields into 
account, a set of masses for the 0++,0"+,l++,l"+,2++ and 2"+ glueball states are obtained and 
are in fairly good agreement with the predictions given by the lattice simulation In addition, some 
new glueball states are predicted. 

PACS numbers:11.10.Qr, 12. 38. Aw, 12.40.Bx, 14.80.Er. 
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o ■ 

O ■ I. INTRODUCTION 

^. 

Oh! -A^s mentioned in the preceding paper, searching for glueballs, nowadays, is an challenging task in particle physics. 
^ ^ Since there are numerous technical difficulties of giving unambiguous identifications of the glueballs in experiment 
^ [1-15]: it is expected that the existence of glueballs and their properties could be precisely predicted by theoretical 
investigations so as to guide the experimental searches. Various methods were proposed in the past to serve such 
^ ' investigations [16-31]. However, the predictions given by different methods are not consistent sometimes and even 
• ^ , contradictory with each other [32-33]. Of these methods, the lattice simulation [26-31] is considered to be faithful. 

■ Nevertheless, even for this method, there still are controversies on the results given by different calculations. Apart 
H ' from the lattice computation, the nonrelativistic potential model [16-18], the relativistic Dirac equation [22] and the 
- - Bethe-Salpeter (B-S) equation [23-25] have recently been applied to evaluate the glueball spectrum. In these methods, 

the interaction between gluons is constructed by two parts: the short-range part which is described mainly by the one- 
gluon exchange interaction and the long-range part which is represented by a phenomenological confining potential. 
In the nonrelativistic potential model, the short-range interaction is simulated by a potential which is derived in the 
approximation of order jt? where v is the gluon velocity and c the light velocity with the assumption that the gluons 
in a glueball move not too fast. In a recent work by using this model [18], with the choice of the lightest glueball 
masses given in the lattice computation [29] as input, the authors obtained a series of two-gluon glueball states with 
masses below ?>GeY . Except for some glueball masses which are in pretty good agreement with the lattice predictions, 
the other calculated masses are apparently different from the lattice results. As was emphasized in Ref.[18], to gain a 
physical solution to the lightest scalar glueball, it is necessary to additionally introduce a phenomenological smearing 
function to replace the (5-function in the attractive contact terms of the potential. Otherwise, the Hamiltonian would 
be unbounded from below. This probably is an unnatural feature of the nonrelativistic potential model. In Ref.[22], 
the calculation of the glueball spectrum was performed by employing the relativistic Dirac equation and showed only 
three theoretical results for the lowest glueball states 0++, 2++ and 3++ some of which are not in so good agreement 
with those given by lattice investigations. In the calculation, a Fermi-Breit potential (the t-channel one-gluon exchange 
potential) was inserted into the Dirac equation with the assumption that the nature and the force between two gluons 
are the same as between two quarks. It seems that this assumption ignores the difference between the potential 
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for quarks and the one for gluons. In addition, it would be mentioned that the Fermi-Breit potential is derived in 
the nonrelativistic approximation of order u^/c^. Therefore, the calculation is not fully relativistic. The relativistic 
calculation of the glueball spectrum was carried out in the framework of B-S equation [24,25]. Owing to the difficulty 
of solving a relativistic equation, only a few states were predicted in these calculations. It is noted that in all the 
previous applications of the B-S equation, the four-dimensional equation was recast in a three-dimensional form in 
the instantaneous approximation in which the retardation effect is completely neglected. Another point we would like 
to note here is that in the aforementioned works, the gluons are all viewed as massive. Each of such gluons in general 
has three degrees of freedom of polarization. Correspondingly, a gluon field should includes three independent spatial 
components: two transverse fields and one longitudinal field in the three-dimensional space. In this sense, we can 
say that the Coulomb gauge as taken in Ref.[25] is inappropriate for the massive gluons because in this gauge the 
longitudinal mode of the field is completely eliminated. Similarly, in Ref.[24], only the transverse gluons are taken 
into account even though the temporal gauge adopted in the work allows existence of the longitudinal gluons. 

In this paper, we intend to investigate the glueball spectrum based on the three-dimensional relativistic equation 
for two-gluon bound states which was derived in our former paper in the angular momentum representation. This 
equation is actually a coupled set of equations satisfied by the four B-S amplitudes for a glueball state: one is related 
to the positive energy states of two gluons, the other three are related to the two gluon states in which there is at least 
one gluon in the negative energy state. In the next section, we will derive from this coupled equations an equivalent 
equation obeyed by only the B-S amplitude of the glueball state for which the two gluons are in the positive energy 
states and give the effective interaction Hamiltonian in the equation a complete form. Since we are unable to compute 
all the terms in the Hamiltonian at present, we are limited ourself to work in a semi-phenomcnological model by which 
the interaction Hamiltonian in the equation is given by the one-gluon exchange kernel plus the phenomenological linear 
confining kernel as was usually done in the previous literature [16-18,23-25]. The new aspects of this paper which 
distinguish from the previous works are: (1) The calculation is fully relativistic and hence includes the contribution 
arising from all the relativistic effects to the glueball masses; (2) The retardation effect of the one-gluon exchange 
interaction is completely taken into account; (3) Apart from the transverse modes of the gluon fields, the contribution 
from the longitudinal mode of the field to the glueball spectrum is appropriately considered; (4) The renormalization 
efi^ect is considered by the effective QCD coupling constant which was derived in the one-loop approximation and 
in a mass-dependent subtraction in our previous work [34]. This coupling constant is not only suitable in the high 
energy domain, but also in the low energy regime; (5) We work in the angular momentum representation. In the this 
representation, the glueball states are easily constructed. In particular, with completing the radial integrals containing 
three and four spherical Bessel functions, the gluon vertices are given explicit and analytical expressions which greatly 
facilitate the numerical task of solving the equation. The theoretical results obtained in this calculation are in quite 
good agreement with those given in the lattice study [29-31], In addition, some new predictions are presented. 

The remainder of this paper is organized as follows. In Section II, we will derive the three-dimensional equation 
satisfied by the gluon positive energy state B-S amplitude from the coupled equations derived in the preceding paper. 
In Section HI, the interaction Hamiltonian obtained in the tree diagram approximation will be discussed and its 
explicit expression will be given. Section IV serves to derive the expression of the linear-wise potential which is used 
to simulate the gluon confinement and incorporated in the glueball equation for numerical calculations. In the last 
section, the calculated results will be presented and some discussions will be made. In Appendices A and B, the 
analytical expressions of three-line and four-line gluon vertices are derived respectively. 

II. THE THREE-DIMENSIONAL EQUATION FOR THE GLUON POSITIVE ENERGY STATE 

AMPLITUDE 

The three-dimensional equations derived in the preceding paper (which will be called paper I later on) are shown 
in the following 



where E is the total energy of a glueball state, uja and uip represent the energies of free gluons 1 and 2 respectively, 
Xafiin) stands for the B-S amplitude describing the glueball state which is defined by 



and K{af3; pa; E) designates the interaction kernel whose closed expression was derived in paper I. In the matrix 
notation, it is of the form 




(2.1) 



(2.2) 
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3 3 3 

= {Y,9iSi - 9iGi393 + 9iGiG-^Gjgj}S-\ (2.3) 

where the matrices in the above expression were clearly defined in paper I. Noticing the definitions of Wa and (see 
Eqs.(4.9) and (5.16) in paper I), 



(fc) ifCa = l 



and 



a+(fc) if^„ = -l ' ^ ' 

where the subscript a on the right hand side (RHS) of Eq.(2.5) is defined without including and hence aa{k) 
and a+(fc) represent the annihilation and creation operators respectively, the equation in Eq.(2.1) may be separately 
written as 

[E„ - w(fci) - (^{k2)]Xa+i3+{n) =J2Kia+(3+;pa;En)xpa{n), 

per 

[En - oj{ki) + uj{k2)]Xa+0- {n) = E K{a+(3-;pa; En)Xpa{n), 

K +a;(fci) -a;(fc2)]Xa-/3+(n) =Y:K{a-p+;pa;En)XpAn), ^^'^^ 

pa 

[En+uj{ki) +w{k2)]Xa-f3-{n) = T,K{a-p-;pa;En)Xp<7{n), 

pa 

where the superscripts ± in and denote ^a)^/3= ±l,the indices p and a still include the indices and ^o- 
respectively and 

Xa+p+{n) = (0+|aaa^|n),Xa+/3-(") = (0+|aaa+|n), 
Xa-/3+(»^) = (0+|a+a/3|n),Xa-/3-(n) = (0+|a+a^|n). 

Following the procedure described in Ref.[35] for fermion systems, the coupled equations in Eq.(2.6) can be reduced 

to an equivalent equation satisfied by the B-S amplitude Xa+p+(n) for the glueball state in which each of gluons is in 
its positive energy state. For later convenience of derivation, we define 

Aab{E)=E-au{ki)-buj{k2), (2.8) 

where the subscript n in En has been suppressed, a,b = ±1, 

(j)++{a(i]E) = Xa+(i+{n),(j)+-{a(3]E) = Xa+p-{n), ,^ gx 

^-+(a/3;i;) =Xa-/3+("),</'— (a/3;£;) =Xa-/3-(") ^ " ' 

and 

Kabcd{al3; pa; E) = K{a''(3''; p'a"; E), (2.10) 

in which the a,f],p,(7 are defined without including the index ^. With the definitions given in Eqs.(2.8)-(2.10), the 
equations in Eq.(2.6) can compactly be written as 

Aab{E)(l>ab{ap;E) = J2J2^'^'>cd{al3; pa; E)<j,,d{pa; E), (2.11) 

cd pa 

where a,b,c,d = ±1. In the product space of momentum ki^ k2 and angular momentum marked by a, the above 
equations may be written in the matrix form 

A++{E)cP++{E) = K++++{E)cP++{E) + J2 K++,a{.E)4>,d{.E), (2.12) 

cd#++ 
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Aab{E)<f>ab{E) = Kat,++{E)cj>++{E) + Kabcd{E)<i>od{E), (2.13) 

where ab ^ ++ and the terms related to have been separated out from the others. Furthermore, In the space 

spanned by 4>ab{E) with ah ^ ++, we use the matrix representation defined as follows 

^{E)=cl>++{E),^{E)={cl>ab{E)}, 
A+{E) = A++{E),A{E) = {A^b{E)}, 

K+{E) = K++++{E),k\e) = {K++,a{E)}, (2.14) 

G{E) = {Kab++{E)/Aab{E)},}, 

G{E) = {K,bcd{E)/A,b{E)}. 
According to these definitions, Eqs.(2.12) and (2.13) may be written in the full matrix form 

A+{E)i;{E) = K+{E)iP{E) + K\E)(t>{E), (2.15) 

(j){E) = G{E)iIj{E) + G{E)(f){E). (2.16) 

Solving the equation (2.16), we obtain 

m = i_^(^) GW(ig)- (2.17) 

Substituting the above expression into Eq.(2.15), we finally arrive at 

A+{E)ij{E) = V{E)i^{E), (2.18) 

where 

V{E) = K+{E) + k\e)^-^G{E), (2.19) 
which is identified itself with the interaction Hamiltonian. Noticing the definition 

3-1— = ^G(")(i;). (2.20) 



71=0 



Eq.(2.19) can be written as 



where 



]/(£;) = ^ !/(")(£;), (2.21) 



n=0 



V(^KE) = K+iE)^ 

v(^\e) = k\e)g{e), 

V^^\E) =k\e)G{E)G{E) 



where 



(2.22) 



According to the definitions in Eq.(2.14), Eqs.(2.21) and (2.22) may be explicitly written as 

V{,aM5;E) = ^l/(")(a/3;7<5;i?), (2.23) 



V(°)(a/3;7<5;£;) = K{a+ (3+]^+5+;E), (2.24) 



E — auj(ki) — buj{k2) 
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_ ■sp ■ip K(a+l3+;p''a'';E)K(p'^a'';ij.''v'^;E)K(tj,''v'^;i+5+;E) 

— 1^ 1^ 1^1^ (B-aw(fci)-6a;(fc2))(B-cw(fei)-dw(fc2)) (2.26) 

067^++ cd!^++ pa- p.v 



where we have used the notation shown in Eqs.(2.8) and (2.10). From Eqs.(2.19)-(2.26), it is clear to see that the 
negative energy states of two gluons act as intermediate states to appear in the effective interaction Hamihonian and 
give contributions to the higher order terms in the Hamiltonian. The equation (2.18) written in an explicit form is 
such that 

[E - uj{ki) - a;(fc2)]V(a/3; E) =Y, V{ap; jS; E)iP{jS; E). (2.27) 

7(5 

This is just the equation satisfied by the gluon positive energy state amphtudes for the glueball states in which 

V(a/3;£)=Xa+/3+(n)- (2.28) 



III. THE HAMILTONIAN GIVEN IN THE LOWEST ORDER APPROXIMATION 



In this section, we plan to discuss the interaction Hamiltonian in the tree diagram approximation of the order of 

here g is the QCD couphng constant. This Hamiltonian can only be given by the term shown in Eq.(2.24) because 
the other terms in the effective Hamiltonian give the contributions which are all higher than . In general, the 
K {a^ (3^ ] ^'^ 5'^ ] E) should be calculated according to the expression denoted in Eq.(2.3) which includes three parts. 
The last part in Eq.(2.3) plays the role of cancelling the B-S reducible diagrams contained in the first two parts and 
gives no contribution of the order g'^. This is because (1) the coefficients 51 (a/3; per A) and gz{a(i] paX) are proportional 
to g and g2{oi(3; parX) is proportional to g"^; (2) the Green's functions Gi{pa\;^5]ti — 12) and Gz{pa\;^5;t\ — 12) 
vanish in the lowest order approximation. As for the first part in Eq.(2.3), in the lowest order approximation, it is 
easy to verify that 

Si{pa\^5) = (0|[: apa^ : aA,a^aa]|0) = 0, 

S2{p(TT\,-i5) ^ (0|[: apa^ar : aA,a^a5]|0) = 0, (3.1) 
53(paA,75) = (0|[: c+c„ : aA,a^a5]|0) = 0, 

where |0) denotes the bare vacuum state. Therefore, we only need to consider the second part in Eq.(2.3). In this 
part, the term related to g2{af3; parX) gives the contribution of order g"^ in the lowest approximation and hence is 
beyond our consideration. For the terms associated with 513(0/3; pa A), the relevant Green's functions vanish in the 
lowest order approximation. Thus, we are only left with terms in the second part of Eq.(2.3) such that 

K°{a+f3+;^+6+;E)=J2Ha+P+;pa;E)S-\pan+S+), (3.2) 

pa 

where 

Aia+(3+;pcr;E) 

= - E E gi{a+p+;^ri>^)G,i{^vX,p,.T;E)g,{^i,yT;pa) (3.3) 

and the indices p,a, ■ ■ ■ should be understood as p^, cr^, • • •. 

Let us first compute the inverse S'~^(pcr; 7+(5"'"). For this purpose, we operate on the both sides of Eq.(3.2) with S 
from the right and get 

A{a+0+;pa;E)=^K''{a+P+;^+6+;E)S{j+6+;pa). (3.4) 

7(5 

It is easy to verify that except for S{^^S~^; p~a~), the S{'-f~^S^; p~a^), S{'j~^S^; p'^a~) and S{'-f~^S~^; p^a"^) are all 
vanishing in the lowest order approximation. As for the 5(7"'" (5+ ; p~a~), we have 

5(7+ (5+; p~a~) = 6jp5sa + S^a^Sp- (3.5) 
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Substituting Eq.(3.5) in Eq.(3.4), we find 

A{a+(3+;p-a-;E) = K°{a+P+;p+(7+;E)+K^{a+l3+;a+p+;E). (3.6) 

Since wo may interchange the indices p and a in Eq.(2.1) or in Eq.(2.27), noticing Xprri'iT-) = Xapin) or ij){pa;E) = 
il){ap;E), we can write from Eq.(3.6) the following relation 

K°{a+f3+;p+a+;E) = ^Aia+p+;p-a-;E), (3.7) 



which means that we may set 



S-\pa;'y+6+) = S-\p-a-;'y+5+) = ^S^pSsa- (3.8) 



Combining Eqs.(3.3) and (3.7), we have 

K°{a+p+;p+CT+;E) 

= -5 E E giia+0+nSX)Giii^S\;fii^T;E)g^{pi^T;p-a-). (3.9) 

■ySX HVT 

In accordance with the definition of gi{ap;'y6X) and g\{p,yT; pa) (see Eqs.(5.7), (5.9), (5.10) and (6.4) in paper I), we 
can write 

51 («+/?+; 7(5A) = E./l(7<5'^)Aa+/3+;TA, 

gi{pvT; p- a-) = E /i(/xJ^A)Ap-^-.;^^, ^^'^^^ 

X 

where 
and 

fi{af3j) = A{af3^) + A(a7/J) + Ai^yap), (3.12) 

here A{a(3j) is the three-line gluon vertex given in the angular momentum representation. Considering the expressions 
in Eq.(3.11) and the fact that only A^+ff- = — A„-^+ = Sa^ are nonvanishing for Aaj3, Eq.(3.10) can be represented 

as 

ffi(a+/3+;7(5A) = /i(7'5a")VA + /i(7'5/?")^a+A- C3 

On inserting Eq.(3.13) into Eq.(3.9), one gets 

K\a+(3+-,p+a+;E) = ^[A{a+ P+ ; p+ a+ ; E) + {a ^ p) .3 

+ {p <r-^ a) + {a <r-^ p,p ^r-^ a)], ^' ' 

where 

A{a+p+;p+a+;E) = ^ h{^da-)GnhSf3+ ; pua' ; E)h{fiiyp+) (3.15) 

and the other terms in Eq.(3.14) can be obtained from the first term by exchanging the indices as shown in Eq.(3.14). 

Now, let us calculate the Green's functions Gii{jS(3'^ ; fiucr' ; E) in the lowest order approximation which arc the 
Fourier transform of the Green functions Gn (7(5/3+; /xi^cr^; ti — ^2)- With the aid of Wick theorem, it can be found 
that only the following Green's function is nonvanishing 

Gn(7+5+/3+;/i-!/-a-;*i-t2) 

= (0|T{: B^{h)Mti) ■■ Mh) ■■ a+(i2)a+(t2) : a+(i2)}|0), ^"^"^^^ 
where cLy(fi) and a+(t2) are the annihilation and creation operators in the interaction picture. Noticing 
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a^(ti) = a^e-''^-*%a+(t2) = a+e^'^<'*^ (3.17) 
and applying the Wick theorem, we find 

Gii{j+6+P+;ii-u-a-:ti-t2) 

= 0{h - t2yr'^'^-'^^'^'^'^^^''-'^\^^,y.h.?>!i. + ^^M^h'^ (3-18) 
By the Fourier transformation and using the famihar integral representation of the step function, we obtain 



(3.19) 



= I dte^^'Gnil+S+P+;n-i.-a-;t) 
Substituting the above expression in Eq.(3.15), we are led to 

j-^ En-u:-,-us-uifs+ie'l '-y I /lU^V; K J (3 20) 

+m-l-p+)]50a + [fi{S-(3-p+) + Mp-6-p+)]Sj, 
+[/i(7-/3-p+) + h{l3-j-p+)]Ssa}. 

Observing the above expression of A(a+/3+; p^a^; E), we see, the first term containing S/^a is an unconnected term. 
It gives the one-loop correction to the gluon propagator whose effect will be included in the QCD effective coupling 
constant by the renormalization procedure. The second term proportional to 6jcr describes the one-gluon exchange 
interaction between two gluons. The third term is the exchanged term for the one-gluon exchange interaction. By 
dropping the first term, we have 



A{a+f3+;p+a+;E) 

li 

E—ojt —Wo- —ojp+ie 



^ ^ [/i(T+<T+a-)+/i(^+r+a-)][/i(T-/3~p+)+/i(/3-T-p+)] ^ (3.21) 



where the summation over r is performed with respect to the gluon intermediate states and the function /i was 
represented in Eq.(3.12) in terms of the function A{af3^) whose explicit expression is derived in Appendix A and 
shown in the following. 



A{a,a2a^) = -|(f )i/«fcfc3 n -^fc,BS*(Z,),,,,T,3,^. 



>^Jl[l'^l'^{kl,k2,k3)T{l^,l'^,mi,^i), 
where 



3 

1 k-fliiflA. _ T, , 

(3.22) 



_ „ ^2ni-i[-i (3.23) 

E '^2(^i+/i2+/i3),ii+'2+'3 n r(ui+i)r(ui-(;+i) ' 

l,/i2, (1*3=0 1=1 " ' " ' 



i n (-l)('*+'*+™*+i)-"[^^l[(2Zi + 1)(2/^ + 1)]5 

1=1 





1'2 




( h k h \ 









I mi m2 ms / 



(3.24) 



r;^3A^ and B^* {h)xiTi are defined respectively in (A. 5) and (A. 15). It is noted that for a given set of angular momenta, 
due to the restriction of S2{fj,^+fj,2+tj.3)j[+i'^+i'^, only a few terms in the series of Eq.(3.23) survive. 
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IV. THE GLUEBALL EQUATION WITH INCLUSION OF THE CONFINEMENT 



As shown in Sec. II, the interaction Hamiltonian in the exact relativistic equation for the glueball states is much 

comphcatcd. In practical investigations, usually, one only considers the lowest order term in the Hamiltonian which 
was explicitly derived in the angular momentum representation in the preceding section. Obviously, in order to get 
reasonable theoretical results, it is necessary to introduce a certain confining potential to phenomenologically simulate 
all the higher order terms in the Hamiltonian [16-18, 23-25]. How to determine the form of confining interaction in 
the relativistic equation ? In this paper, as was similarly done for the quark-antiquark system [36], we introduce the 
confining Hamiltonian operator in such a manner 



^ / d''xid}'x2r"" r'^'A^ i^i) ■ x^i)y(|x^i - a^2|)3^( x^2) • A^ilt^), (4.1) 



H, 

__ 1 
~ 2 

where 'A^{'x') stands for the gluon field operator and V^d"^! — "^2!) the confining potential which will be taken to be 
the linear one [37] 

V(|^i-^2|)=7l^i-^2|, (4.2) 

here the parameter 7 designates the strength of the confining potential. When the ghion field operators in Eq.(4.1) 
are replaced by their expansions in terms of the multipole fields (see the expansion given in Eq.(4.14) in paper I), 
Eq.(4.1) will be represented as 

Hc= ^ Vc{aia3;a2a4) : (4.3) 

aia2a3a4 

where 

Fc (a 1 as; 0:2 0:4) 

Here the symbols Ai and Pi were defined in Appendix A. This is just the wanted confining potential written in 
the angular momentum representation which will be inserted into the relativistic equation. By using the Fourier 
transformation 



(4.5) 



J {2ny -ct^ 
the expansion for the plane wave function 

e^-^ = 4n^i^Mpr)Yi*^{p)YUx) (4.6) 

Im 

and the expression for the scalar multipole field 

Af^iklt) = ^J^kji{kr)Yim{x), (4.7) 

we can write 

f°° 1 

y(|^i-^2|) = -87r75^ / dq-^Af^{q-^i)AC{q^2). (4.8) 

Im ^ 

Substitution of this expression into Eq.(4.4) leads to 

ye(aia3;a2a4) = -47r7r^7='^'=^ / dq-D^,^,rD*^^^^,, (4.9) 

where 
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Dca^r = J d^Xi^^lil^i) ■^^l{l>^i)Af{^i), ^^^^^ 

here r = (q, Z, to). 

Completely analogous to the calculations described in Appendix A, it is easy to derive the following expression 



Im i—1 

^Ji[i'2iiki,k2, q)Jiy'j{k3, ki, q) 
xT{liJ'i,rrii,r]i,l,m,l)T{lj,l'j,mj,r]j,l,m,-l), 



(4.11) 



where i = 1, 2 , j = 3, 4 



T{liJ'^,mi,r]i,l,■m,r]) 

= n(-i)(''+'^+™'+i)'^''^[^^T^i(-i)"*«'"[^^^^i ^'^''^'^^ 

xT{li,l[,riimi,l,rim) 



with 



T{li,l'^,mi,l,m) = / dn{x)y^i^i>^.rn^{x) • T^j^i^m^ (x)y;m(S) 

l[ 1'2 l \ f h h I X\ h k I 



X 1 I -\ 

y I mi m2 m } \1'2 I'l 1 



(4.13) 



and Ji.i',i{ki, kj,q) is defined as the same as given in Eq.(3.23). With the introduction of the above confining potential, 
the total interaction Hamiltonian in Eq.(2.27) is now taken to be 

V{ap; jS) = Vg{ap; jS) + Vdap; jS), (4.14) 

where 

Vg{al3;^6)=K%al3;-rS), (4.15) 

which was formulated in Eqs.(3.9)-(3.24) and the Vc{a(3;^d) was given in Eqs.(4.11)-(4.13). 

Now we turn to discuss the wave function ip{ap) in Eq.(2.27) which was defined in Eq.(2.28). In the lowest order 
approximation, the two-gluon bound states can be written in the form 

where 

/a/3 = "y|'^ciC2C';{mii2TO2/Aili,A2i2('''l'^2), (4-17) 

in which (5ciC2 represents the color singlet, Cj^^^_^i^^^ is the C-G coupling coefficient, a = (ci, Ai, Zi, mi, fci, = +1); 
(3 = {c2,X2,l2,m2,k2,^i3 — +1), J, M are the total angular momentum and its third component of a glueball and w 
denotes the spatial parity and charge conjugation parity. With the introduction of the cluster coordinates 

I^ = ti + t2, t = i(ti-t2), (4.18) 

where 1^ and It are the total momentum and relative momentum respectively, we see, in the center of mass system 
(i? = 0), Eq.(4.17) reads 

fSp = -^6,,c.Cf^Z,i^^MXMi.ikm, - k2), (4.19) 
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where k = \ k \ = ki = k2. Substituting Eq.(4.16) into Eq.(2.2), we find 



= = ^^c,c.QX;.m.5A,t.A.i.(fc)'5(fci - k2), (4.20) 

where 

9Xlh,\2l2ik) — fXlll,X2l2 + i~^)^f\2l2,Xlll' (4-21) 

in which h = l\ + I2 — J- Evidently, if Ai = A2,/i = I2 and h = odd, we have g'y^ll Aab^^) ~ ^' "^^^^ gives a new 
selection rule for the glueball states. Upon substituting Eq.(4.20) into Eq.(2.27) and noticing 



a = (c, Xi,h,mi,k,+l),P = (c, X2,l2,m2,k,+1), 
p = (c', A3, 13, ma, q, +1), a = (c', A4, h, rui, q, +1) 



(4.22) 



and 

" " dk, (4.23) 



a cXlm ' 

we finally arrive at 



{E - 2oj)g(^ll ^ I {k) 

= mfdqV{X^hk; Mhh Xskq; XdiT, E)gif^^^j^{q), (4-24) 

where oj = \Jk'^ -V p?, E is the mass of a glueball state given in the center of mass frame and 

V{Xilik; A2/2A;; As/ag; XiUq; E) 

^("/3;P^;i^)C,X;.m.^?;lm3Um., (4-25) 

mim3m2m4Cc'M 

here V{af3-, pa: E) was given in Eq.(4.14). When the explicit expression of V{af3-, pa: E) is substituted in Eq.(4.25), 
one can see that the summation over 7711,7712,7713,7714 and M is easily carried out by utilizing the well-known formula 
for the angular momentum coupling and the summation over the color indices c and c' can be completed by noticing 
jabcjabc _ 24. We think, it is unnecessary to show here the result given by these summations. 

The equation in Eq.(4.24) is the eigenvalue equation used to calculate the glueball spectrum. In the calculation, the 
QCD coupling constant g contained in the part of Hamiltonian V^(a/3; 7^) is replaced by the running one which was 
derived in Ref.[34] recently in the one- loop approximation and in a mass-dependent momentum space subtraction. 
The coupling constant used in this calculation is of the form 



where as{X) = g^{X)/A-iT, 



G(A) = 11 In A - '^Nf[2 + ^/Stt - ^ + (1 + (4.27) 



in which Nf is the number of quark flavors. 



r,(A) = ^/xP^ln ^ ^ ^ (4.28) 

and A = J^^P^ here ■p is chosen to be the transfer momentum of the exchanged gluon which may simply be taken 

V QCD 

as p = k — q (ov simplicity and Aqcd is the QCD scale parameter. The running coupling constant shown above is 
applicable not only in the high energy domain, but also in the low energy regime. Particularly, in the large momentum 
limit, it immediately goes over to the result obtained previously in the minimal subtraction scheme. 



10 



V. NUMERICAL RESULTS AND DISCUSSIONS 



In this section, we first show the theoretical glueball masses calculated from the equation given in Eq.(4.24) and then 
make some discussions. Our calculation is performed by using the standard program of Mathematica which allows 
us to compute the effective Hamiltonian in Eq.(4.24) analytically. In this paper, we confine ourself to investigate 
the low- lying glueball states including the 0++,0 '",1++,! '",2++ and 2 ground and lower excited states whose 
masses are less than A.OGeV. Some of these states have been investigated before in various models. We also examine 
the effects of the longitudinal mode of the multipole fields and the different sets of free parameters on the glueball 
masses. In our calculation, the theoretical parameters are adjusted so as to be able to compare our results to those 
presented recently by the lattice simulation [29]. The parameters taken are: the gluon mass = 0A2GeV which 
is comparable with yu = (0.5 ± 0.2)GeV taken previously in the nonperturbative continuum studies[38], the scale 
parameter Aqcd = OABGeV and the strength of the confining potential 7 = O.lSGeV^ , which satisfies the relation 
/Lt ~ Aqcd and 7 ~ Ag^^^ which make the parameters essentially depend on a single dimensional quantity. Moreover, 
the value of 7 = O.lSGeV^ is consistent with that of the string tension in lattice simulations. The coupling constant 
= 0.3 and quark flavor Nf = 3. The calculated masses of glueball states are displayed in table I. In the table, the 
case I and the case II respectively denote the results obtained with and without considering the contribution arising 
from the longitudinal mode of the multipole fields which appears in the intermediate states of the matrix elements 
of the interaction Hamiltonian. In the last column of the table, we quote the results shown in Ref.[29] which were 
calculated by the lattice simulation. 



Table I. The mass spectrum of two-gluon glueballs. 



Glueball states (J^*^) 


Mass(GcV) 
Case I Case II 


Mass(MeV) 
Lattice results 


0++ 


1.73 2.18 1730(50)(80) 
2.66 3.59 2670(180)(130) 
3.59 


0-+ 


2.60 2.30 2590(40) (130) 
3.65 3.78 3640(60)(180) 


1++ 


2.73 2.42 
3.45 3.51 


1-+ 


2.67 2.59 
2.87 3.01 


2++ 


2.43 2.43 2400(25)(120) 


2-+ 


3.32 2.26 3100(30)(150) 



As seen from Eq.(4.24), each glueball state is not only assigned by its spin J and parity tt, but also related to the 
mode marked by (Ai, li)(A2, ?2)- In this paper, we take low-lying modes to perform the calculation. For the scalar 
glueballs of quantum numbers J^'-' = 0++ and the tensor ones 2++, according to the angular momentum and the 
parities of the multipole fields, we take the mixture of the modes (TElTEl) and (TMlTMl). For the glueball states 
^ and 2 ^, the modes arc taken to be (TElTMl) for every glueball, as was similarly done in the investigation 
within the bag model [19,20]. This means that these glueballs are mainly constructed by the gluons with transverse 
polarization. But, this does not imply no contribution of the longitudinally polarized gluons to these glueballs. The 
longitudinal gluons may, as virtual particles, appear in the intermediate states in the effective interaction Hamiltonian. 
It is emphasized here that for the transverse mode of gluons, as mentioned in Appendix A, the mode h = is not 
permitted. This mode can only exist for the longitudinal gluons. Different from the case of massless gluons, the 
longitudinal mode of massive gluons is possible to take part in formation of some glueballs. For example, the glueball 
states 1 can be formed not only by a combination of modes (MlE'l) and (SILO) which gives the states with masses 
as listed in the table, the modes (MlEl) and (LILO) can also form the glueball states with masses 3.2306^ and 
3.82GeV respectively. For the states l++,we only take the mode (LOMl) in our calculation because according to the 
B-S amplitude constructed in Sec.IV (see Eq.(4.21)), the modes (ElEl), (MlMl) are forbidden. 

In order to determine the parameter dependence and errors in our calculation, we take three sets of parameters in 
which we set /x = Aqcd and 7 (in unit GeV^)= Aqqj-). The results of case I in table I with different parameters (the 
parameters and Nf remain unchanged) are presented in Table II. We find that the masses of glueballs increase 
gradually when the set of parameters increase. This indicates that our numerical calculation is stable and the results 
are reliable. 

Table II. The mass spectrum with different parameters. 
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Mass(GcV) 








0++ 


i.DD 


1 7/) 
i. /4 


l.oi 






9 fiS 


9 




O.OZ 


O.OU 


O.DD 


n-+ 




9 fin 

^ .uu 


9 fiS 




3.57 


3.66 


3.73 


1++ 


2.64 


2.73 


2.82 




3.37 


3.46 


3.54 


1-+ 


2.59 


2.69 


2.78 




2.78 


2.88 


2.96 


2++ 


2.35 


2.45 


2.53 


2-+ 


3.25 


3.34 


3.41 



The calculated results show that the gluon mass give an appreciable effect on the glueball masses. In our calculation, 

the mass of gluon should be aroimd 0.45 GcV. This fact indicates the reasonability of the QCD with massive gluons 
which is chosen to be the starting point in our calcidation. For this kind of QCD, it is necessary to take the longitudinal 
mode of gluons into account. As shown in Table I. when the longitudinal mode is considered in the intermediate states, 
the theoretical masses for the states mentioned above would be greatly improved. Otherwise, there would occur a 
considerable discrepancy between the results given by this paper and the lattice calculation. In addition, as illustrated 
before, the longitudinal mode allows us to investigate more glueball states which possibly exist in the world. 

In comparison with the previous theoretical glueball masses obtained from the B-S equation and the Dirac equation, 
our results give more support to the lattice predictions [29-31] which are believed to be more reliable because the lattice 
calculation is based on the QCD first principle and essentially nonperturbative. Within the statistical errors existing 
in the lattice calculations, our results shown in the first column of Table I can be considered to be consistent with 
lattice predictions for the low-lying glueballs with masses less than 3.5GeV, especially, for the lowest scalar glueball 
state 0++ with mass about 1.7GeV and the tensor glueball state 2++ with mass about 2AGeV. The achievement 
of the better consistence is obviously attributed to the fact that our calculation is fully rclativistic and is able to 
include the contributions arising from the retardation effect and longitudinal mode of the gluon field which could 
not be considered in the previous investigations [23-25]. Now let us analyze our results in some more detail. It is 
mentioned that the lowest scalar glueball state 0+"'" and the tensor one 2"'"+ have been investigated in many models 
and the theoretical masses are almost the same even though there are a little difference between different calculations 
[29-31] . For example, the mass of the lowest state 0++ was given by 1754 ± 65 ± 86MeV in a recent lattice calculation 
[31] which is different from that given in Ref.[29]. It is expected that these states may be identified with the pure 
glueball states and searched out first in future experiments. Aside from the two states mentioned above, the first 
radial excited state 0++ with mass 2.66GeV should correspond to the state 2670(180)(130)Mey given in the lattice 
simulation even though whether the latter is a pure glueball or not is still in question [32,33]. The next radial excited 
state 0"*"^ with mass 3.59GeV is a new prediction given in this paper which was not predicted in the lattice simulation 
and the other calculations. As for the pseudoscalar states ^ and pseudotensor state 2 ^ are all comparable with 
the corresponding states presented in the lattice calculation. But, the mass of the state 2 ^ is little higher than the 
lattice one. In addition, we note that the states 1""*" and 1++ were not predicted in the lattice simulation, but the 
states 1 with masses 2.67GeV and the state 1+"*" with mass 2.73GeV are compatible with the recent calculation 
by the nonrelativistic potential model [18]. 

In conclusion, it is emphasized that different from some previous investigations, the calculation in this paper is based 
on the rigorous three-dimensional relativistic equation satisfied by the two-gluon glueball states which is derived from 
the QCD with massive gluons and represented in the angular momentum representation. Especially, the interaction 
Hamiltonian in the equation is given a complete expression which provides a firm basis for further study. In this paper, 
even though we work in the relativistic potential model with introducing phenomenologically a confining potential, the 
new consideration of the retardation effect and the longitudinal mode of the gluon fields allows us to get the improved 
theoretical results which are well consistent with the lattice predictions. The only uncertainty in our calculation arises 
from the introduction of confining potential. Certainly, if a sophisticated confining potential could be found from the 
exact interaction Hamiltonian derived in this and former papers, it would be anticipated that a relativistic calculation 
may give more accurate theoretical predictions. 
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VII. APPENDIX A: THE EXPRESSION OF GLUON THREE-LINE VERTEX 



This appendix is used to derive the explicit expression of the gluon three-Hne vertex in the angular momentum 
representation. As shown in Eqs.(4.7), (4.29) and (4.30) in paper I, the gluon three-line vertex in the interaction 
Hamiltonian is 



ij3 = -^pbc j d^x(A^ X I^) • (V X A^). 
It is represented in the angular momentum representation as follows 

Hg= ^ A{aia2a3) : aa^aa^aa^ ■, 

where 

Aia.a^as) = -^r''^ J d^'xil^ x 1^) ■ (V x 



(A.1) 



(A.2) 



(A.3) 



Here we have set a, = {Xi,f3i) in which Aj = TE,TM, L mark the transverse electric, transverse magnetic and 
longitudinal modes of the multipole fields and /3j = (Zj, m,, fcj, ^j). From now on, we use the symbols lirui to represent 
the total angular momentum. For later convenience, the relations between the multipole fields which were mentioned 
in paper I are represented in the matrix form 



V X 



Aim [kX j 



1 
fc I 1 




/ it'^^ I- — >^ 



(A.4) 



When we define 




(A.4) can be rewritten as 



Thus, (A.3) reads 



With the definition 



Vx Ai;^{k^)=kTx,x,Ai;^{k^). 



(A.5) 



(A.6) 



(A.7) 



(A.8) 



where I' = I + t and r = 0, ±1, the relations shown in Eqs.(3.2)-(3.4) in paper I for the multipole fields can also be 
written in the matrix form 



—^TM 



1 



V 



/ -^'C\k-t)\ 



(A.9) 



If we define {TE, TM, L) = (-1, 0, 1) and 

B{1) = 



10 



(A.10) 



13 



where I ^ and 




B{0) =10 



-i 

which means only the longitudinal mode survives when / = 0, then (A. 9) can concisely be written as 

After inserting (A. 12) into (A. 7) and noticing the definition given in Eq.(4.16) in paper I, we have 

A(aia2a3) 

where we have defined 



yiVm^kl^)-\ tu'^(klt) ifC = -l 



and 



^^'>^--\ B*{h)^, ife--l ■ 
In light of the the expression in (A. 8), the integral over in (A. 13) can be represented as 



where 



and 



Ji[i'j'^ (fci, fc2, fcs) = J drr^jv^ (fcir)j;^ {k2r)3v^ {k^r). 



On substituting (A. 16) in (A. 13), we just give the formula denoted in Eq.(3.22). 
First, let us calculate the T{li, rrii, ^i) in the case of = 1. In this case, we set 

^{li,l'^,mi) = j dn{x)^ i^ii^^^{x) x T*i,;^„,(x)] • 'fi^v^^^ix). 

By using the following formulas [39,40]: 

^lVm{x)=Y.C\Z,,,^Yi.^,{x)l^,, 
m'q 

6 g = (""I)'' 6 _g, 6 g ■ ^ q' — ^qq'i 

we find 



nm^nzqiq^qs 

X / dQ.{x)Yv^ni {x)Yi>^n2 {x)Yv^m {x). 
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Employing the familiar formula for the above integral and the definition and property of 3-j and 9-j symbols for the 
angular momentum couplings [39,40], one can get 



_ ,r(2ii+i)(2f2+i)(2i3+i)(2;'i+i)(2;^+i)(2;f,+i)ii 



2tt 



l[ 1'2 I'Af h h h\ \ u u u I (A.22) 

In the case of ^ = —1, it is easy to find 

tl^ix) = {-iy+''+^+'tu--m{S). (A.23) 



If we define 



V^" _ / '^ll'mix) if ^ = 1 



the vector spherical harmonics ^ii'm{x) and Y^;;/„(x) can be represented in an unified form 

tl^ix) = (A.25) 
Thus, according to (A.22)- (A.25), we have 

T{li,l'i,mi,r]i) = f dn{x)[f]'i,^^^{x) x i^l','^„,(x)] • '^Z'.mA'^) 

= n(-i)('.+':+™^+i)-[^ir(^„;^,r,,mO, ^^-^^^ 

i=l 

where T{li,l'^,ri^mi), as expressed in Eq.(3.24), is directly written out from (A.22) with replacing m,i by rjimi. 

Let us turn to compute the integral in (A. 18) following the method proposed by one of the authors in this paper and 
his coworker in their early publications [41]. As we know, there is a momentum conservation in the gluon three-line 
vertex: ki + /c2 + hi = which gives a certain restriction on the magnitudes of the three momenta. In fact, from 
ki = { k2+ /cs)^ = fc| + fc| + 2fc2^3 cos 012, it is seen that when cos0i2 = ±1, we have ki = k2 + ks and fci = |A:2 — fel. 
This implies that only when the conditions fci + ^2 ^ ks,k2 + k^ > fci and ki + k^ ^ k2 or 

ki + k2-kz> 0, fc2 + fcs - fci > 0, fci + fcs - fc2 > (A.27) 

are simultaneously satisfied, the momentum conservation holds; whereas, when ki > k2 -\- k^, the momentum conser- 
vation is violated. In addition, adding any two inequalities in (A.27) together, we find ki > 0, therefore, each ki varies 
from zero to infinity. In later derivations, the following relations are useful 



hf\x) = ji{x) + ini{x), hf\x) = ji{x) - ini{x), ,^ 

ji{-x) = i-iyjiix), h^\-x) = {-iyhP{x), 



where ji{x) is the spherical Bessel function, ni{x) the spherical Neumann function, h'l^\x) and h^'^\x) arc the first 
class spherical Hankel function and the second class one respectively. The asymptotic behaviors of these functions 
are as follows. When x ^ 0, 

^,(^\ , 2'V.x' „ /^N ^ -(2^-1)!! 

Jl[X) > (2i+l)!> ni[x) > ^i + i , 

uWr^] , -'(20! 1,(2)/ N , i(2i-l)!! y^-^^) 



and when 



First, we prove that the integral Ji'^i;^i',^{ki, k^, fca) vanishes in the case of fci > fc2 + fca. In this case, considering the 

analytical property of the functions ji{x) and h^i'\x) as shown in (A. 29) and (A. 30), the following integral along the 
contour C on the upper half complex plane of r as depicted in Fig.l is zero 



15 



£ drr''h'il\hr)ji,{k2r)ji,{ksr) = 0. (A.31) 

The contour C can be divided into four parts, C = Co + (— oo, 0~)+Ci + (0"'", +00). Clearly, the integral along the large 
half circle Ci vanishes when |r| tends to infinity. Thus, noticing the relations in (A. 28), Zj > and I1+I2 + I3 = even 
which is implied by the first 3-j symbol in (A. 22), one can get from (A.31) 

JiW,{ki,k2,h) = -l j drr^hll\kir)ji,{k2r)ji,{hr). (A.32) 

J Co 

Substituting the series expansions 



^ly'^'^i - 2 r(M+i)r(/x+i+|)' 

^ °° (■_j^-)i+^+i(^-)2M-!-i (A. 33) 

ni{kr) = — E r(^+i)r(^'-(+i) 



into the right hand side of (A.32), it is easy to find that the integral along the circle Co around the origin also vanishes 
when \r\ goes to zero. Thus, we reach the following result 

Ji[i',i!,{ki,k2,ks) = 0. (A.34) 

Next, we compute the integral under the conditions shown in (A. 27). In view of these conditions and the asymptotic 
behaviors of h\^\x) and /ip'(a;) shown in (A. 30), the function /(r) defined by 



fir) = h\l\k,r)h\'Jik2r)h\'J{ksr) + hf^\k,r)h\f {k2r)h^^J {k,r) 
+h^^\k^r)h\l\k2r)h[l\k,r) + h'il\k,r)h'il^ik2r)^^^^ 



is analytical on the upper half complex plane of r except for at the origin. Therefore, we have 

drr^fir) = 0, (A.36) 



X 



where the contour C is still represented in Fig.l. Due to the conditions in (A. 27), the integral along C\ still vanishes. 
Thus, from the above integral, we get 

JiW,{kuk2, fca) = drr^fir) + drr'f{r)} 



In accordance with (A. 28), the function /(r) can be written as 

(A.38) 



m 

= 4jji {kir)ji.^ {k2r)ji3 (ksr) + 2ini, {kir)ji.^ {k2r)ji^ {k^r) 
+2iji^{kir)ni^{k2r)ji^{k'ir) + 2iji^{kir)ji^{k2r)ni^{k3r) 
+2ini^ {kir)ni^ {k2r)ni^ (ksr). 



Inserting this expression into (A. 37) and using the series representation in (A. 33), it can be found that except for the 
last term, the other terms in (A.38) all give no contribution to the integral. Therefore, we have 

Jl[i'^i'^{ki,k2,k3) = J^^drr'^ni,{kir)ni^(k2r)ni.Jksr) 

_ iEl V n i-^r-^'H^)'^'-'--' r , 2(^i+^,+M3)-;i-i.-i3-i (A-39) 
- 32 2^ 11 ^{|J.i+l)^{^n-l^+^) Jco 

Setting r = pe^^ and noticing 2(/xi + 112 + Ms) — h — h — h = even and /° d9e^^"^^ = —TT5m,o here m is an integer, we 
finally obtain the expression as shown in Eq.(3.23). 
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VIII. APPENDIX B: THE EXPRESSION OF GLUON FOUR-LINE VERTEX 



In this appendix we would like to derive the explicit expression of the gluon four-line vertex in the angular momentum 

representation for completeness although the vertex gives no contribution to the equation (2.27) in the lowest order 
approximation due to 52 (/cutA, 7^) = as shown in Eq.(3.1). The four-line vertex in the interaction Hamiltonian at 
t = which was described in Eqs.(4.7), (4.29) and (4.31) in paper I may be rewritten as 

ij4 = aLjabcjcde J ^^^.(ji . A^){^ . A^) 

where the second equality is obtained by substituting the expansion of gluon fields in terms of the multipole fields 
into the first equality and the coefficient function is of the form 

By making use of the representation in (A.8)-(A.12), (A. 14) and (A. 15) and noticing the definition given in Eq.(4.16) 
in paper I, the B{aia20i3a4) can be represented as 

V 1 2 J 4; -^j J ./^iMM) ^ ' (B.3) 



where 



and 



Ji'^i'^i'^U {ki, k2, ks, ki) = j drr^jr^ {kir)jr^ {k2r)ji'^ {ksr)jv^ {kit) (B.4) 

f{k, I'i, rrii, rji) = j rfO(J)T^J'^Vimi (^) • ^w>2 (^)^W>3 (^) • ^w>4 (^)' (^-5) 

here the notation in (A. 24) has been used. Inserting (A. 20) and (A. 25) into (B.5) and employing the familiar formulas 
for the integrals of spherical harmonics and for the angular momentum coupling [39,40], it is not diflicult to get 

f{li,l'i,mi,r}i) 

= — X;im(-1)''"^'^"^'=''^'''''^"(2^ + l)f]1^_^(-l)(''+'i+"'' + l)'^'°[ "T'" l 



X 



^4^2 1 J V Vi'mi VsiTT'S rn J \ 7727712 7747724 -m 



To compute the integral in (B.4), we first examine the conditions satisfied by the magnitudes of the momenta. From 
the momentum conservation fci + fc2 + fcs + ^4 = 0, we have 

fc? = (^2 + fcs + klf = kl + kl + kl ,g ^^ 

+ 2fc2fc3 cos 023 + 2fc2fc4 COS (?24 + "i-k^k^ COS 634. 

From the above equality, we may find the maximum and minimum of ki by setting Oij = or tt. There are 
only four ways which permit us to take the values of 9ij = or tt. These are (1) 623 = O24 = O34 = 0; (2) 
623 = ^24 = TT, 634 = 0; (3) 023 = ^34 = TT, 624 = 0; (4) 024 = 034 = TT, 023 = 0. Correspondingly, we get from (B.7) the 
equalities ki = k2 + k3 + k4, ki = \k2 — k3 — k4 \ , ki = \k3 — k2 — k4\ and ki = 1^4 — ^2 — ^sl ■ From these equalities 
we find the following restriction conditions which are consistent with the momentum conservation 



k2 + k3 + k4- ki > 0, fci + fcs + fc4 - ^2 ^ 0, 

ki + k2 + k4 - fca ^ 0, /ci + ^2 + ^3 — ^4 ^ 0, ,„ ^\ 

ki + k2 - k3 - k4 ^ 0, ki + k3 - k2 - k4 ^ 0, \ ■ ) 

fci + ^4 - ^2 - fcs > 0, fci + ^2 + fcs + ^4 > 0. 



17 



By adding some two of the above inequalities, one may sec ki ^ 0. And similar to the proof described in the preceding 
appendix, it can be proved that the integral in (B.4) is merely nonvanishing provided that the conditions in (B.8) are 
respected. According to the above conditions, it is obvious that the function defined below is analytical on upper half 
complex plane of r 

F{r) 

= h['^ {hr)h\'^ {hr)h\'^ {hr)h[f [k^r) + ft« (fcir)/i« {k2r)h^i'J {k3r)h\'^ (^4^ 

(fcir)/i|f ik2r)h'il'^ ik3r)h['J (kir) + /if ik,r)h\'J [k^M''} {ksr)h[^ (k^r) (B.9) 

+hll^ {k,r)h['J {k2r)hf^ {k,T)hf^ ik,r) + /i« (fcir)/z|f (fc2r)/z« {ksr)h<i^^ {k.r) 
{k,r)hf^ {k2r)hf^ {k,T)hf^ {k,T + {k,r)hf^ {k,T)hfJ {k,r)hfj {k,r) 

where the have been replaced by li for convenience. Therefore, based on the Cauchy theorem, the integral along 
the contour C as depicted in Fig.l is zero, 

/ drr^F{r) = (B.IO) 
Jc 

Prom this equation, noticing h+l2 + h + li = even as implied by the first two 3-j symbols in (B.6), we obtain 

JiwM^iMMM) = ^{j,^drr^Fir)+j'^^drr'F{r)} 

In view of the relations in (A. 28), the function F{r) can be represented as 



F{r) 

= 8ji^ {kir)ji^ {k2r)jis {k^r)^^ (kir) + 2ji^ {kir)ji^ {k2r)ji^ {kir)ni^ {kAr) 

+2jii {kir)ji^ {k2r)ni^ {k?,r)ji^ [kir) + 2iji^ {kir)ni.^ {k2r)ii, {ks.r)i^ (k^r) 
+6mii {kir)ji^ {k2r)jh {k3r)ju (kir) + 2ini^ {kir)ji2 {k2r)ni^ {k3r)n^ {kiv) 
+2ini^ {kir)ni^ {k2r)ji^ {k3r)ni^ {k^r) + 2ini^ {kir)ni^ {k2r)ni^ {k3r)ju {k4r) 
-2iji^ {kir)ni^ {k2r)ni^ {k3r)n^ {k4). 



(B.12) 



Upon inserting the above expression into (B.ll), using the scries representation in (A. 33) and considering the 
relations among the angular momenta which are implied by the 3-j and 9-j symbols in (B.6), it is easily verified that 
the first five terms in (B.12) give vanishing contributions to the integral. The nonvanishing contributions given by 
the last four terms can be calculated by the same procedure as described in the last part of Appendix A. The result is 

Ji[i'^ii,uiki,k2,k3,k4) 

o , CXD 

— Z!_('-'l5('l+'2+'3+U) V A . . „ „ „ „ 

— ley ) 2^ '^2{tii+H2+H3+H4.),l'.2+l'3+l'i-l'i /T^1o^ 

A'l,M2,A'3./J4=0 \p.L6) 

(1-2^.1) yr (k./k^f^^-'j 

^ ^ r(^i+i)r(/*,+;.+|) . 11 fe,r(/*,+i)r(/*,-;,+i) ' 
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X. FIGURE CAPTION 



Fig.l. The contour for the integrals containing three and four spherical Bessel functions. 
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